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Abstract 

A factorial-moment analysis with real (integer and non-integer) phase space par- 
tition is applied to vr^'^p and K"''p collisions at 250 GeV/c. Clear evidence is 
shown for self-affine rather than self-similar power-law scaling in multiparticle 
production. The three-dimensional self-affine second-order scaling exponent is 
determined to be 0.061±0.010. 



" Supported by the Government of the Republic of Armenia, contractnumber 94-496 
Permanent address: Department of Physics, Jingzhou Teacher's College, Hubei 434100, China 
Onderzoeksdirecteur NFWO, Belgium 

Supported by the Polish State Committee for Scientific Research 



1 



1 Introduction 

Classical non-Abelian theories show non-linear behavior and the non-linear development 
of a classical non-Abelian system can be determined numerically PI. Although extension to 
a quantum theory as QCD is far from trivial, non-linearity is inherent also to QCD shower 

development J^] 

The first experimental evidence for non-linear behavior in high-energy multiparticle pro- 
duction came from multiplicity fluctuations in a JACEE event recorded in 1983 I^'^l. In this 
event, the total multiplicity is about one thousand and the fluctuations in small rapidity bins 
arc 2 times average. In 1987, NA22 t^l found an event in which the fluctuation in a small 
rapidity bin is as high as 60 times average. 

To be able to decide whether these fluctuations are dynamical, i.e. larger than expected 
from Poisson noise, Bialas and Peschanski [^1 suggested to use factorial moments (FM), defined 
as 

1 ^ {rimirim - 1) ■ • • (n^ - g -H)) 
^'^'^ = Mi^^ ■ « 

In the above equation, M is the partition number of a fixed phase space region A under 
consideration, 6 = A/M is the size of a sub-cell, rim is the number of (charged) particles 
falling into the mth sub-cell. The authors show, if the power-law scaling 

FgiS) oc S-'l'^ (2) 

holds when 6^0, then dynamic self-similar fluctuations are present in multiparticle produc- 
tion. 

In the following, the behavior of Fq in ever smaller phase-space cells S was studied in 

almost all high-energy experiments on lepton-lepton, lepton-nucleon, hadron-hadron and 
nucleus-nucleus collisions and approximate scaling was established (for recent reviews see 
[''^1). It has been pointed out 1^1, however, that the anisotropy of phase space has to be 
taken into account and that the multiparticle final state in a high energy hadron-hadron 
collision may be 'sclf-affine' ^^^^ rather than self-similar. 

In this note we apply a self-affine analysis to the NA22 data on vr+p and K+p collisions 
at 250 GeV/c, where, in addition, the phase-space partition M of (1) is generalized to real 
(i.e. also non-integer) values. This method reduces the x^/NDF values by a factor of 2 
with respect to the limitation to integer values and grants clear evidence for the presence of 
dynamical sclf-affine fluctuations in these collisions. The three-dimensional self-affine scaling 
exponent is determined to be (j)2^ = 0.061 it 0.010. 

2 The method 

A self-affine transformation in the three phase-space variables denoted Pa^Pb^Pc is defined as 
^Pa ^Pa/^a-, ^Pb ^ ^Pa/\-, ^Pc ^ ^Pc/ ^c, with shrinking ratios A^, and Ac, respectively. 
The anisotropy (self-affinity) of a dynamical fluctuation can then be characterized by the 
so-called roughness or Hurst exponents 1^°] 

^^ii = {i,j = a,b or a,c or 6,c), (3) 

with Xi<Xj, 0<Hij<l. (4) 
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These exponents can be obtained '^^1 from the experimental 6 dependence observed in the 
second-order factorial moments in the three corresponding variables. If self-affine fluctuations 
of multiplicity exist in multiparticle production, exact scaling, i.e. a straight line in InFg 
versus InM, should be observed if and only if the Aj are allowed to differ from one another. 

The above prediction was checked on our data [12] and on 400 GeV/c pp data.l^^] ^j^g 
Hurst exponents for longitudinal-transverse directions were determined to be 0.474ib0.056 in 
the {y,(p) plane and 0.477ib0.057 in the {y,PT) plane for 250 GeV/c 7r"'"p and K"'"p collisions 
and 0.74ib0.07 in the (77, (p) plane for 400 GeV/c pp collision. A self-affine higher-dimensional 
analysi^ was performed in both cases and the results are confirmative. 

However, the original analysis was limited by the small number of combinations of (inte- 
ger) partitions Mi allowed by condition (3). Therefore, the method (and necessary correction 
procedure) has recently been generalized to non- integer partitions M.t^^l In this letter the 
new method and correction procedure are applied to obtain a precise answer on the question 
of self-affinity in our data. 

2.1 Non-integer FM analysis [^''1 

To be definite, let us consider a one-dimensional analysis in rapidity y. In the ideal case, the 
factorial moments Fq{5y) depend on the bin width 6y = ^yjM^ but not on the position of 
the bin on the rapidity axis. If that is the case, the result of averaging over all M bins as in 
(1) is equal to that of averaging over bins with N < M. This means that ideally one has 
F2(M) = F2{N,M), where 

F2{N, M) = 1 f: Mn^- 1)) ^ <M,6y = Ay/M) . (5) 

This equation can be used as the definition of FM for any real (integer or non-integer) value 
of M [1^] , provided that the number N of bins used for averaging is taken to be 

N = M-a, (0<a<l). (6) 

However, even in the central region the rapidity distribution is not flat. The shape of this 
distribution influences the scaling behavior of the FM. Therefore, the cumulant variable 

was introduced, [1^1 which has a flat distribution by deflnition. An additional correction factor 
has to be introduced again for the FM analysis with non-integer partition. 

To see this, let A denote the phase-space region in consideration, 6m the mth bin, pi(yi) 
and P2{y 1,112) the one- and two-particle distribution functions, respectively. Then we have 

{um) = f pi(2/)(iy = — /" dy2 [ ^2/1^2(2/1, 2/2); (8a) 
(f^m(ram - 1)) = / dy2 j dyip2{yi,y2)- (86) 



^In the case of NA27, due to lack of momentum measurement, the self-affine analysis was 2D instead of 
3D. 
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After transforming to the cumulant variable, (n^) becomes a constant, independent of 
m. However, comparing the two above equations, it can be seen that due to the difference in 
the integration region over y2, {nm{nm — 1)) is in general not constant even though (n^) is. 
This was experimentally verified on our data. 

Note that in the definition of -Fg in (1) a horizontal average is taken. When the partition 
number M is an integer, the horizontal average is over the full region A. The variation of 
{nm{nm — I)) / {nrn)'^ IS thus smeared out, and no correction is needed. On the contrary, when 
M is non-integer, the horizontal average is performed over only part of the region and the 
influence of the variation of {nm{nm — l))/(^m)^ becomes essential. 

2.2 Correction factor for the — l))/{nmy distribution 

From the definition of F2{N,M) in (5) it can be seen that only N bins arc included in the 
horizontal average when M is non-integer (M = N + a, 0<a<l). Consequently, only a 
fraction r = N/M of the full region A is taken into account. We minimize the influence of 
this by introducing a correction factor R{r) and define t^^l 

In order to extract R{r) = R{N/M) from the experimental data, let us go back to integer 
M and calculate F2 averaging only over N of the M bins (A'' < M). The result is, in general, 
a function of both N and M. The correction matrix 

N 

C{N,M) = ^^ , iV = l,2,...,M. (10) 

m=l 

is shown in Fig. 1. as a function of N/M for M = 3,4, . . . ,40. Points for different M lie in 
a narrow band (mind the scale), so that the correction factors Ry{r), i?p^(r) and R(^{r) can 
be obtained from an interpolation of C{N, AI). 

It turns out, however, that the result is very sensitive to the interpolation function. An 
inappropriate choice of the interpolation function, even inside the narrow C{N, M) band, will 
be either insufficient in eliminating the "sawteeth" observed in the uncorrected lni<2 versus 
InM plot or it will over-correct. It is found l^^l that when the "sawteeth" lie above the smooth 
curve of integer M, as in the case of y and pT, the upper boundary of the C(N, M) band has 
to be used for the interpolation. When the "sawteeth" reach from above the smooth curve 
of integer M to below, as in the case of ip, the middle of the C{N, M) band has to be used. 
The interpolation functions used for the three cases are shown as dotted lines in Fig. 1. 

Having obtained the correction factor for one-dimensional FM's with non-integer M, the 
correction of higher-dimensional FM's can be obtained from a straight-forward generalisation 
[-^^1 paying special attention to the overlap regions (Mi — Ni)(M2 — N2) and (Mi — A^i)(M2 — 
N2){M^ — N3) in the two- and three-dimensional analysis, respectively. 

3 Results 

The details of the EHS spectrometer can be found in '^^^ , those of the trigger and data analysis 
in [^^1. The acceptance criteria and data samples are those already used in I^^l. 
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3.1 One-dimensional analysis and Hurst parameters 

In Fig. 2 are shown the results of a one-dimensional analysis of lnF2 versus InM. The first 
column reproduces the results obtained earlier for integer M l^^l. The second column gives 
the results for real M as defined in (5). It can be seen that the points for non-integer M 
depart in a "sawtooth pattern" from the curve determined by the points with integer AI, 
especially in the cases of y and pT- After correction by Ry, i?^^ and i?^ obtained from the 
interpolation of the correction matrix C {N, M) above, the results (shown in the third column 
of Fig. 2) become smooth. 



3.2 Two-dimensional analysis 

The plots for 2-dimensional self-affine FM's are presented in Fig. 3 for -f/prp^, = 1 and Hypr^ = 
Hy^ = 0.475. The same 1-dimensional correction factors Ry, Rp^ and R,p are used together 
with the geometrical factors taking care of the overlap regions l^^l . The corrected results for 
{y,PT) and (y, (^) are satisfactory. For {pt,v) -^2 is improved considerably with respect to 
the original real (integer and non-integer) My plot shown in the middle column of Fig. 2. 



3.3 Three-dimensional analysis 

The final 3-D results are presented in Fig. 4. The left figure is the one with integer M 
obtained before ^^^^ The right figure is the corrected result of the self-affine analysis with 
real M together with a linear fit 

lnF2 = A + BlnMy . (11) 

In order to minimize the influence of momentum conservation [^^1, the fit starts from My = 2. 
For the same reason, the non-integer M points are given for My > 2, only. 

It has to be noted that integer M values could be used in our previous analysis only 
because the Hurst exponents had values close to 0.5 and 1.0, respectively. This was rather 
accidental, however, and a generalization to real (integer and non-interger) partition t^^l is 
necessary to be able to evaluate other experiments (e.g.^^^l). A comparison of integer-M and 
real-M results in table 1 shows that the slope B itself is unchanged but the error tends to 
decrease. Furthermore, it can be seen from table 1 that the ratios x^/NDF are decreased as 
compared to those obtained from integer M and at the same time the number of degrees of 
freedom (NDF) increases largely. 

Note that the slope B of (11) is related to 02 of (2) via 



MsD = MyMp^M^ = My 



as 

02 = i?/fl + -^ + -^') . (12) 

From the B value for weighted real M in table 1, we get as the final result for the scaling 
exponent 

= 0.061 ± 0.010 . (13) 
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4 Conclusions 

In this paper we presented the results of a self-affine analysis of experimental data on factorial 
moments generalized to real (integer and non-integer) partition M. Correction factors were 
introduced to minimize the influence of the variation of {n,n{n„i — 1)). The corrected results 
for non-integer M lie on smooth lines interpolating between the integer-M points obtained 
before. 

The results of a three-dimensional self-affine analysis are well fitted by a power law, giving 
a more general and more precise check of the self-affine power-law scaling of the data than 
the previous analysis based on integer partition only. The observed behavior is consistent 
with the fact that the longitudinal direction is privileged over the transverse directions in 
hadron-hadron collisions. It would be important to use this analysis in the study of QCD 
parton-shower development in e'^e" collisions. 
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Table 1. The parameter values obtained from a fit of the 3D-data by (11) 



Method 


A 


B 


XVNDF 


Without bin-size correlation 








weighted integer-M 


-0.04±0.03 


0.32±0.03 


7/4 


weighted real-M 


-0.01±0.02 


0.32±0.02 


19/20 


Without bin-size correlation 








unweighted integer-M 


-0.08±0.02 


0.33±0.03 


12/4 


unweighted real-M 


-0.05±0.02 


0.33±0.02 


30/20 


With bin-size correlation 








unweighted intcgcr-M 


-0.08±0.02 


0.34±0.02 


14/4 


unweighted real-M 


-0.05±0.02 


0.34±0.02 


51/20 
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Figure captions 

Fig. 1 The correction matrix C{N, M) as function of N/M for M = 3 - 40 (for the curves 
see text). 

Fig. 2 The one-dimensional plots of lni<2 versus InM. The first column are the previous 
results for integer M. The second column are the results of real M as defined in (5). The 
third column are the results after correction. For easier comparison, the first six integer-M 
points are indicated as full circles in columns two and three. 

Fig. 3 The same as Fig. 2 for two-dimensional self-affine FM as a function of InM^. 
Fig. 4 The three-dimensional plots of lni*2 versus InM. The left figure gives the previous 
results for integer M. The right figure gives the corrected results of the self-aiRne analysis 
with real M together with a linear fit. 
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Fig. 3 



